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ABSTRACT 



The generation of gravity waves in a frontal zone is investigated 
with a numerical model. The solutions show that gravity waves can be 
generated in a frontal zone by the application of a heating function 
which varies sinusoidaly with time, and that such waves propagate both 
ahead of and behind the frontal zone. The amplitude of the waves are 
directly proportional to the magnitude of the forcing function and to 
the period of sinusoidal oscillation. Even with large heating rates, no 
hydraulic jumps are formed. 
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I. ■ INTRODUCTION 



Meteorologists involved in numerical weather prediction have for some 
time been filtering gravity waves out of their prediction models. This 
has been done because high frequency inertial-gravity oscillations could 
be generated which would give unrealistic pressure changes, obscuring the 
meteorological forecasts. Gravity waves may be meteorologically signifi- 
cant however in explaining irregularities in the ionosphere, and in some 
mesoscale phenomenon in the troposphere. Hines (1960) has suggested that 
internal atmospheric gravity waves may account for many of the irregulari- 
ties observed in the lower ionosphere and that these waves may have their 
origin in the large energy regions of the lower troposphere. Gossard, 
Jensen and Richter (1971) have studied gravity waves of small amplitude 
with a high-resolution FM/CW radar sounder. Breeding (1972) investigated 
the creation of unstable regions in the atmosphere by internal gravity 
waves. Gossard (1962) has listed the following three mechanisms for 
generating internal waves in the troposphere: 

”1. Internal waves can be generated as standing waves in the lee of 
topographic features. 

”2. Gravity waves are often caused by a traveling boundary between 
two types of air of different density. This type of tropospheric wave 
train is short, lasting less than two hours. It is a fair-weather 
phenomenon, rarely associated with storm fronts. 

**3. Oscillations which last for many hours are generally associated 
with tropospheric storms and synoptic scale features. These oscillations 
of long duration are usually generated by storms and frontal disturb nces 
but the presence of a stable layer in the lower atmosphere is requi^^cid." 
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Tepper (1952) postulated that gravity waves may be generated in the 
vicinity of a moving frontal zone if the zone moves into an area where 
a stable layer exists. As the waves move out they steepen into hydraulic 
jumps. He felt that large vertical motions caused by the jumps could 
lead to severe weather such as is found near squall lines. Kurth (1971) 
studied hydraulic jumps in a rotating, continuously stratified atmosphere. 
Numerical simulation of squall line formation was carried out by Sasaki 
(1959) and Ogura and Charney (1962). One of the objectives of the present 
study was to determine if such effects could occur near a properly simulated 
frontal zone. 

Since the advent of satellite pictures and the technique of producing 
film loops showing ”weather-in-motion” it has become evident that high 
speed gravity waves move out from the face of frontal zones (Fujita and 
Bohan, 1967). 

In this study, the non-geostrophic frontogenesis model of Williams 
(1972a) is modified in an attempt to generate gravity waves in a frontal 
zone. 

In section II, the Boussinesq equations are given and the modeling 
relations are introduced. In the model the time-dependent quantities 
are functions of y and z only. The basic wind deformation field is 
independent of time and height. The numerical procedures, boundary 
conditions and initial conditions are given in section III. Numerical 
results are presented and discussed in section IV. Conclusions and 
recommendations for further investigations are given in section V. 
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II.. BASIC EQUATIONS 



In this experiment the Boussinesq equations are employed and the 
domain is bounded by two rigid horizontal planes. The compressibility 
of the atmosphere, which is neglected in the Boussinesq approximation, 
should not be of qualitative importance since the density scale height 
in the atmosphere is much larger than the thickness of typical frontal 
zones. The replacement of the tropopause by a rigid horizontal surface 
will induce large errors in this region and prevents gravity waves from 
propagating out of the troposphere, but the resulting errors near the 
lower boundary should be small. Basic equation development follows the 
work of Williams (1972a). 

The hydrostatic Boussinesq equations can be written in the following 
form when the earth’s rotation and horizontal eddy viscosity are included, 

^ V + V- (w) (wV) + V4> + /■ k V = V, (II- 1) 

V- (eV) t = + Q, (II. 2 ) 

V- V + = o , (II.3) 

= 3©/@, , (II-4) 

where is a constant reference potential temperature; p^ a constant 
reference pressure, k = R/ (Cp) ; T(po/p)^ - 0 ^ , the departure of 
the potential temperature from 0^, (j) = Cp 0,(p/p 6)^4- gz, the pressure 
function, A^ ^ V , the horizontal diffusion of momentum, 0, the 

horizontal diffusion of heat-, and Q the heating function. 
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An exact steady solution to this set of equations is given by 



V = D (xi-yj) 
w = 0 

^ = — D^(x^+y^) / 2 J -f Dxy 

0 



, (II. 5) 



where D is a constant, Cartesian coordinates have been introduced 
and f is taken to be constant. The horizontal velocity is a field of 
pure deformation and this deformation is given by 2D, 

If departures from the fields (11,5) are independent of x, then they 
will remain independent of x. Thus, the dependent variables are subdivided 
as follows: 

V =j^D x+ u(y,z,t)j i+ -Dy+v(y,z,t)j j J 

w = w(y,z,t) I 

^(y,z,t) I . (II. 6) 

+ r (y, z, t) ‘ 

If these expressions are substituted into Eqs, (II, 1)- (II, 4) , they 



become: 

^ +. ^(^ 4 ) , y/u a. VIH - -/Zr - fl /TT 7^ 

4^ » (II. 10) 

ay ^ 

, (II. 11) 

©. 

where V = -Dy, The boundary conditions are 

. w = 0, z = 0, H , (11.12) 



where H is the distance between the horizontal plates. 
If the vertical average of a quantity is defined as 

<(( ^ 
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and the hydrostatic equation (II. 11) is integrated with respect to z 



and the vertical mean is removed, then 






(11.13) 



The y equation of motion may now be written in terms of ^ • 

If the vertical average of (II. 10) and the boundary conditions (11.12) 
are used, then 



i, (-> ' 



O 



This equation states that the total disturbance y mass flux is independent 
of y. If the other variables have proper symmetry it follows that ^ 

must be an odd function of y which leads to 

^ • ( 11 . 14 ) 

If the vertical average of (11.8) is taken, and (11.12) and (11.14) 
are used , then 

h<y-vy ^ 4./'<'u> = o 

• (11.15) 

This result is subtracted from (II. 8) which yields 



+ <u» 



(11.16) 

Eqs. (II. 7), (II. 9), (II. 10), (11.13) and (11.16) form a complete set 
which can be solved by a pure marching process. 
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III. NUMERICAL PROCEDURE 



The arrangement of variables and the finite difference approximations 
are the same as those used by Williams (1967). In order to close the 
problem computational boundaries must be introduced in y. Since the 
disturbance velocity should die out at a sufficient distance from the 
axis of dilatation, then 

+ Y, 2, t) = 0 . (III.l) 

However, there is appreciable inflow across these computational boundaries 
since V(+ Y) = + DY. The quantities u and ^ which are advected across 
the boundaries must be specified independent of the interior values if 
computational stability is to be maintained (Platzman, 1954). Thus, the 
following boundary conditions are used: 

u[+ (Y+iQy/2,z,t] = u[+(Y+Ay/2), z, o] 
e[± (Y+Ay/2,z,t] = ©[+(Y+Ay/2), z,0] 



1 



(III. 2) 



The computational boundaries y = + Y are placed between grid points 
so that the above conditions are actually applied at y=+ (Y+Ay/2). 
These conditions are replaced by time specified functions in some 
experiments. 

These boundary conditions developed a weak nonlinear instability 
near the boundaries but this was removed by introducing a forward step 
every 72 time steps. 

The initial temperature field, from Williams (1972a), is given by 
® ~ “H/i) — ^ (x/vr )<?rc^cr/? ]^(III.3) 



33 



-/ — 

were <=^ ~ H ^ quantity 

which is constant, is the initial stability and A is one-half the total 
horizontal temperature variation. The initial x component of the velocity 
is given by 

The initial afield is obtained from the quasigeostrophic equation. 
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IV. NUMERICAL SOLUTIONS 



The finite-difference forms of Eqs. (II. 7), (II. 9), (II. 10), (11.13) 
and (11.16) are solved by a pure marching process. Experiments are 
conducted by varying boundary conditions and heating functions. The 
Williams model is allowed to run for four days by which time a frontal 
zone has stabilized as described by Williams (1972b). At four days, 
boundary conditions are changed or heating is introduced and the model 
runs for an additional four days. The parameters which describe the 
various experiments are given in Table I. 

The procedure for initializing the heating function Q is as follows: 

At the end of four days, the vertical velocity field is scanned at 

z=2.5 km (level three) to determine the area of maximum vertical velocity, 

which corresponds to the leading edge of the frontal zone. A region R 

is defined which is the region where the heating function is applied. 

This region varies in width and vertical extent depending on the 

experiment. The function is applied as follows: 

Q = C g(t) G(y,z) in region R, 

Q = 0 elsewhere. 

Here C is the maximum heating added and g(t) takes the value 1.0, or 
the following sinusoidal variation 
g(t) = s in (cot) 

where cO is the frequency. The function G(y,z) takes either the value 
1.0 for block heating or the value determined by the sine squared function 



N . 2 ( 'TT [ y - ] ) 

G(y,z) = s in { t "; , ' ' : V 

I width of heating zone f 

In later . experiments , a cooling function is applied to the upper Ic 3ls 
of the model to simulate evaporation. This cooling is only activated 
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TABLE I. Variation of parameters 
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when there is heating and it is zero otherwise. The magnitude of the 



cooling is equal to 20% of the heating below. 

Run number one is a control run with no variation of parameters. 
Maximum vertical velocity and the approximate position of the frontal 
zone is shown in Fig. la. Departure potential temperature, which is 



Fig. lb. The divergent part of the wind at four y grid points is 
shown in Fig. 2a and Fig. 2b. Fig. 3a and Fig. 3b show the variation 
in the departure potential temperature in the control run and the same 
quantity in run number three where a stable layer is introduced at the 
beginning. In both cases, all waves are virtually damped out at the end 
of eight days. 

Run number two tests Tepper’s hypothesis that gravity waves will be 
generated if a frontal zone moves into an area where a stable layer exists. 
A stable layer is introduced at the end of four days ahead of the frontal 



was set equal to 8®K km . No gravity waves or jump zones were noted. 

Run number three is the same as the control run except for the 
introduction of a stable layer at both the +Y boundaries at the beginning 
of the program. As in run number two, the layer is inserted between 



these two levels. By the end of four days this stable layer had 
propagated throughout the model between levels z=1.5 km and z=2.5 km. 
The presence of a stable layer throughout the model did not appear to 
generate gravity waves or create jump zones. 




is the initial stability, is shown in 



zone at the -Y boundary between z=1.5 km and z=2.5 km where 

oc. ,^-l . 




z=1.5 km and z=2.5 km and 




is set equal to 8°K km 



-1 



between 
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O 00 




Fig. la. Vertical velocity w maxima and the position of the 
front ( ). Run number 1. 




Fig. lb. Departure potential temperature at level z=1.5 km. 
Run number 1. 
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(l-”* “) * 




Fig. 2a. Divergent part of the wind v as a function of time at four 
levels; z=.5 km (1), z=2.5 km (3), z=5.5 km (6) and z=8.5 km (9) 

Run number 1. 
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Fig, 2b. Divergent part of the wind v as a function of time at four 
levels; z=.5 km (1), z=2.5 km (3), z=5.5 km (6) and z=8.5 km (9). 

Run number 1. 
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Fig. 3a. Departure potential temperature ^ ’ as a function of time a 
four levels: z=.5 km (1), z=2.5 km (3), z=5.5 km (6) and z=8.5 km (9) 
Runs number 1 and 3. 
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In run number four, a heating function of 10®K day is introduced in 
block form and is 120 km wide, with vertical extent from z=.5 km to 
z=4.5 km, and it does not vary with time. This roughly approximated the 
release of latent heat which would be typical of convective activity in 
a frontal zone. As can be seen in Fig. 4a and Fig. 4b an increase in 
vertical velocity and a definite wave pattern in the departure potential 
temperature ^ ’ ahead of the frontal zone was noted. The gravity waves 
appeared to be damped as the distance from the frontal zone increased. 

An additional result of this experiment was that the heating in the warm 
air greatly strengthened the temperature gradient in the frontal zone. 
This coincides with the work of Eliassen (1959) concerning the formation 
of fronts in the atmosphere. 

In run number five the heating function took the shape of a sine 
squared curve with maximum heating occurring in the center of a 300 km 
wide zone with vertical extent from z=.5 km to z=4.5 km. This gave a 
smoother heating function in the frontal zone. The results were similar 
to run number four except that the vertical velocities and gravity x^7ave 
amplitudes were reduced slightly, as might be expected. 

In run number six the heating function was allowed to vary sinu- 
soidaly with time which approximates the build-up and decay of convective 
activity in the vicinity of a frontal zone. The frequency for a gravity 
wave with a wave length equal to the scale of the heating zone gave a 
period of 9.9 hours. As a result of this variation, vertical velocity 
maxima appeared at Y;i^+1300 km which can be seen in Fig. 5a. Other 
bands of vertical velocity relative maxima appear but only the maximum 
values are shown in this and subsequent figures. Departure potential 
temperature fluctuations were also noted at Y;:^+1300 km as can be ^n 
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Fig, Aa. Vertical velocity w maxima and the position of the 
front ( ), Run number 4. 




Fig. 4b, Departure potential temperature Q ’ at level 
z=1.5 km. Run number 4, 
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in Fig. 5b. More significantly, the gravity waves generated in the 
frontal zone propagated in both directions from the front and extend to 
the boundaries of the model. This can be seen in Fib. 6 which shows 
definite wave propagation in the departure potential temperature field 

at various locations on the Y-^axis* 

In run number seven the heating function was reduced to 3 K day 
in order to determine if the patterns found in run number six would still 
be evident. Gravity wave generation and propagation still occurred 
but the magnitude of the waves was greatly reduced. Vertical velocity 
values also decreased and the maximum values returned to the vicinity 
of the frontal zone. This also was the first run where the heating 
function was raised and extended from z=2.5 km to z=6.5 km. This slight 
modification was made to account for the fact that convective activity 
in cumulus clouds does not extend to the surface. 

In run number eight, the heating function was increased to 30®K day ^ 
in order to increase the nonlinear effects in the gravity waves. Vertical 
velocity maxima remained at Y^+1300 km but the magnitude of the vertical 
velocities were greatly increased. Also the first evidence of the 
formation of standing waves was observed. Well defined bands of vertical 
motion maxima were also noted. 

The heating function was increased to 50®K day ^ in run number nine 
and as can be seen in Fig. 7a and Fig. 7b, regions of strong vertical 
motion were noted throughout the model as well as corresponding variations 
in potential temperature. Standing long waves appear to form with a 
wave length of approximately 360 km. 
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Fig. 5a. Vertical velocity w maxima and the position of the 
front ( ). Run number 6. 




Fig. 5b. Departure potential temperature & ' at level 
z=1.5 km. Run number 6. 
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Fig-. 6. Departure potential temperature as a function of time ; 

four levels: z=.5 km (1), z=2,5 km (3), z=5.5 km (6) and z=8.5 Ion (S . 
Runs number 6, 12 (mod) and 12. 
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Fig. 7a. Vertical velocity w maxima and the position of the 
front ( ). Run number 9. 




Fig. 7b. Departure potential temperature ^ ’ at level z=1.5 kra. 
Run number 9. 
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In order to determine what effect a reduction of the horizontal 
eddy viscosity would have on gravity waves, in run number ten Am and 

were reduced to 7,000m^ sec“^ from 20,000 m2 sec~l. This reduction 
appeared to increase vertical velocities and led to standing waves of 
approximately 130 km, as can be seen in Fig. 8a and Fig. 8b. 

In runs eleven through thirteen, the heating function was held at 
30 ®K day-1 with A„j and A^ held at 20,000m2 sec"^ , and the stable layer 
was altered. In run number eleven the stable layer was removed entirely 
and the areas of maximum vertical velocities shifted from 1300 km 

to + 1200 km. In run number twelve a strong stable layer with 
equal to 12®K km”^ was inserted between z=1.5 km and 
z=2.5 km. The maximum vertical velocities moved almost to the extremities 
of the model, away from the frontal zone. Run number thirteeen was a 
repeat of run number two where Tepper’s assumptions were tested, but with 
the inclusion of the heating function. Fig. 9a and Fig. 9b show much 
the same results as previous experiments with the exception that behind 
the frontal zone, where no stable layer exists, the maximum w was moved 
closer to the frontal zone. An interesting result of run number thirteen 
can be seen in Fig. 10 at Y=-840km and Y=-1140 km where a long wave 
developed at z=8.5 km (level nine). This agrees with Alberty and 
Van Sickle (1969) who found evidence of a mesoscale wave between a 
frontal zone and a severe storm. 

Run number fourteen was another control run which includes all the 
changes made to the original model. It offers an opportunity to make 
detailed comparisons with run number one. Table I contains details of 
parameter changes made between run number one and run number fourteen. 
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Fig. 8a. Vertical velocity w maxima and the position of the 
front ( ). Run number 10. 




Fig. 8b. Departure potential temperature at level z=1.5 km. 

Run number 10. 
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Fig. 9a. Vertical velocity w maxima and the position of the 
front ( ). Run number 13. 




Run number 13. 
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Fig. 10, Departure potential temperature ^ * as a function of tin it 
four levels z=,5 km (1), z=2.5 km (3), z=5,5 km (6) and z=8,5 km (^. . 
Run number 13, 
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Fig, 11a and Fig. 11b show the now coramon profile of w and ^ ’ distribu- 
tions. Fig. 12a and Fig. 12b show the divergent part of the wind v as 
a function of time for four positions on the Y-axis ahead of the frontal 
zone. A phase velocity was obtained from the propagation of the phase 
lines of the velocity field. In this case the phase speed was approxi- 
mately 25nsec”l, This can be compared to Alberty and Van Sickles’ 
mesoscale wave velocity of 14m sec”^ from' observed data. Fig. 13 shows 
an interesting phenomenon which was present in all the cases where 
gravity waves propagate away from the frontal zone. As the waves at 
level three (2.5 km) and level six (5.5 km) propagate, the 0 ’ traces at 
these two levels diverge. At the point in time where they first begin 
to diverge, the 0’ trace at level nine (8.5 km) decreases while the 0 ' 
trace at level one (0.5 km) increases slightly. When the & ' traces at 
level three and level six begin to converge again, a long wave appears 
to be generated at level nine, with a marked increase in The 

traces at Y=-1140 km seem to suggest that this divergent/ convergent 
pattern repeats itself. No conclusions will be drawn concerning this 
phenomenon at this time, but it should be an interesting area for future 
theory and study. 

In runs number 15 and 16 the frequency with which the heating 
function varies with time is altered. In run number 15 a frequency of 
2'rT/52,000 is used corresponding to a maximum every 14.5 hours. This 
led to an unrealistic slope in the frontal zone and an outward displace- 
ment of maximum w, as can be seen in Fig. 14a. The departure potential 
temepratures exhibited wide fluctuations as shown in Fig. 14b and 
Fig. 15. In run number 16 a frequency of 2^/20 yOOO is used which ed 
to a wave period of 5.5 hours. The results of this experiment Wv^ . 
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Fig, 11a. Vertical velocity w maxima and the position of the 
front ( ). Run number 14. * 




Fig. 11b. Departure potential temperature^* at level z=1.5 km. 
Run number 14. 
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Fig. 12a. Divergent part of the wind v as a function of time at four 
levels: z=0.5 km (1), z=2.5 km (3), z=5.5 km (6) and z=8.5 km (9). 

Run number 14. 
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Fig, 12b. Divergent part of the wind v as a function of time at four 
levels: z=0.5 km (1), z=2.5 km (3), z=5.5 km (6), and z=8.5 km (9). 
Run number 14. 
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Fig. 13. Departure potential temperature ^ * as a function of tim 
at four levels: z=0.5 km (1), z==2,5 km (3), z=5.5 km (6) and z = 8.5 Rm 
(9). Run number 14. 
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Fig. 14a. Vertical velocity w maxima and the position of the 
front ( )• Run number 15. 




Fig. 14b. Departure potential temperature 0 ’ at level z=1.5 km. 
Run number 15. 
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Fig. 15. Departure potential temperature 0 ’ as a function of time it 
four levels: z=0.5 km (1), z=2.5 km (3), z=5,5 km (6) and z=8.5 km J) . 
Run number 15, 
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quite similar to the results of run number seven in that gravity waves 
were generated and did propagate but, even with C=30°K day”^ the waves 
generated were of very small amplitude. The maximum vertical motion was 
located in the frontal zone and the departure potential temperature 
profile was similar to that of run number one. 

It was expected that if the nonlinear effects were important, that 
the scale of the waves would become smaller as the waves moved out from 
the source region. In order to further examine these effects, the 
boundaries were moved from 1800 km to 2700 km in runs 17-20, These 
experiments in fact do show a decrease in scale as the waves propagate 
out. If this decrease comes from nonlinear effects then it should be 
more pronounced in cases with larger heating. If run 19 with C=3®K 
day“l is compared to run 20 with C=100° day“l, it can be seen from 
Table I that the maximum w for run 20 is less than 33 times the raaximum 
w for run 19. If a nonlinear scale reduction occurred in run 20 then 
the maximum w in that experiment would be larger than the maximum w in 
run 19, times the ratio of the heating rates. This seems to indicate 
that the observed scale may not be caused by the nonlinear effects. 

Rather the observed scale reduction probably comes from the contraction 
in the basic deformation field. 

A comparison was also made with the results obtained by Kurth (1971) 
who found that in a rotating, continuously stratified atmosphere, 
hydraulic jumps will form if a certain relationship exists between the 
Froude number and the Rossby number. The Froude and Rossby numbers 
were determined for runs seven, nine, fourteen, and twenty, using the 
Froude number F=V/[D(g /dz) 1/2 ), and the Rossby number R^= V/fL, 

with the result that all four cases fell on the jump side of the 
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critical curve F=AR^^, where A is a constant in the range of 6. 5-7. 5* 
Even though jumps should have occurred, it is evident from the previous 
discussion that they did not occur. The conclusion must be that gravity 
waves, by themselves, do not cause hydraulic jumps under the conditions 
of this study. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The generation of gravity waves in a frontal zone has been investi- 
gated numerically. It was determined that the presence of a stable 
layer ahead of a frontal zone was not sufficient to generate gravity 
waves or hydraulic jumps. It was found that the addition of a heating 
function on the face of the front strengthens the frontal zone and 
generates gravity waves, but that these waves do not propagate away 
from the frontal zone. It was found that if the heating function 
applied to the frontal face is allowed to vary sinus oidaly with time, 
that gravity waves will be generated at the front and will propagate 
in both directions from the frontal zone and induce vertical motion 
maxima both ahead of and behind the front. It is believed that such 
motions could cause convective activity and condensation on the less 
stable warm air side. The amplitude of the gravity waves and the 
magnitude of the vertical motion maxima varied directly with the heating 
function and the period of oscillation. A stable layer need not be 
present for gravity waves to generate and propagate, although the 
existence and strength of the stable layer appears to affect the 
distance the vertical velocity maxima are displaced from the front. 

The vertical velocity maxima fluctuated with the scale of the model 
but these fluctuations appear to be caused by the contraction in the 
basic deformation field. The amplitude of the gravity waves generated 
did not appear to be large enough to induce nonlinear effects which 
might lead to the formation of hydraulic jumps or squall lines. Squall 
lines do not appear to be formed by a simple Tepper mechanism but ^ ay 
still be caused by topographic effects. 
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In future studies using Williams model, a smaller grid scale should 
be used in order to obtain the fine detail which may be missed using 
a course grid. The model would be improved if moisture were parameter- 
ized or if it were possible to fix the forcing function based on 
observational data. The latter is not yet possible due to inability 
to measure the heating distribution. The model should be modified 
by replacing the rigid plate upper boundary with a boundary which would 
allow gravity waves to propagate out of the troposphere. It is also 
suggested that cooling be included both above and below the heating 
region. 

Possibilities for further study using ojDservational data would be 
to determine if a correlation exists between stable layer strength and 
displacement of severe storms from a frontal zone, and also whether 
divergence and convergence of potential temperatures curves between 
two intermedieate levels in the troposphere leads to long wave formation 
near the tropopause. 
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